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SUMMARY 

A rapid approximate method of analysis was developed for Both com- 
pressible and incompressible, nonviscous flow through radial- or mixed- 
flow centrifugal compressors with arbitrary hub and shroud contours and 
with arbitrary blade shape. The method of analysis is used to deter- 
mine approximately the velocities everywhere along the blade surfaces, 
but no information concerning the variation in velocity across the 
passage between blades is given. 

In eight numerical examples for two-dimensional flow, covering a 
fairly wide range of flow rate, impeller -tip speed, number of blades , 
and blade curvature, the velocity distribution along the blade surfaces 
was obtained by the approximate method of analysis and compared with the 
velocities obtained by relaxation methods. In all cases the agreement 
between the approximate solutions and the relaxation solutions was 
satisfactory except at the impeller tip where the velocities obtained 
by the approximate method did not, in general, become equal on both 
surfaces of the blade as required by the Joukowski condition. 


INTRODUCTION 

In impellers of centrifugal compressors, part of the viscous losses 
and the phenomena of surge and choke are related to the velocity dis- 
tribution on the blade surfaces. Viscous losses in impellers are 
associated with the boundary layer along the flow surfaces. The growth 
of this boundary layer depends on the velocity variation along the flow 
surfaces just outside of the boundary layer. In particular, if the 
velocity decelerates rapidly along the blade surfaces, the boundary 
layer may separate causing large mixing losses. Also, if the velocity 
at any point along the blade surface is sufficiently greater than the 
local speed of sound, shock losses will result. The choke phenomenon 
occurs when the average velocity between blades is sonic. This average 
sonic velocity is characterized by local supersonic velocities along 
portions of the suction surface of the blade. One possible cause of 
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surge ill centrifugal impellers is the formation of relative eddies on 
the pressure surface of the blade (reference l) . These eddies are 
characterized by negative velocities, opposed to the general flow 
direction, along portions of the pressure surface. In order to analyze 
the performance of centrifugal impellers it is therefore necessary to 
determine the velocity distribution on impeller blades. 

Several methods of analysis that can be used to determine the 
variation in velocity along blades with finite spacing have been 
developed for two- dim ensional incompressible flow (references 2 to 5, 
for example) and compressible flow (references 1 and 6). All these 
methods require considerable labor and therefore are not convenient 
tools for analyzing the performance of an arbitrary impeller design. 

In this report a rapid approximate method developed at the NACA 
Lewis laboratory is presented for both compressible and incompressible, 
nonviscous, two-dimensional flow between blades with finite spacing in 
radial- or mixed-flow centrifugal compressors with arbitrary hub and 
shroud contours and with arbitrary blade shape. The method of analysis 
can be used in connection with an axial -symmetry solution to determine 
the velocities everywhere along the blade surfaces, but no information 
concerning the variation in velocity across the passage between blades 
is given. 

Other approximate methods that are less rapid than the proposed 
method for computing the velocity distribution on blade surfaces in 
impellers of centrifugal pumps and compressors are given in refer- 
ences 6 to 9. In the sections SIMPLIFIED ANALYSIS of refer- 
ences 6 and 7 approximate methods are developed for computing the 
theoretical distribution of velocity across the passage along normals 
to the blade surfaces. The methods are limited to straight or 
logarithmic-spiral blade shapes on radial or conic surfaces of revolu- 
tion and do not apply, because of assumptions, in regions near the 
impeller tip and the impeller inlet. In reference 8 methods are 
developed for computing the distribution of velocity across the pas- 
sage between blades in the circumferential direction for incompressible 
flow with arbitrary blade shapes and with arbitrary hub and shroud 
contours. The methods do not apply, because of assumptions, in regions 
near the impeller tip and the impeller inlet. In reference 9 an 
approx ima te method is developed for computing the theoretical velocity 
distribution everywhere within the impeller. In this method the cor- 
rections required for compressibility and for blade unloading at the 
tip are somewhat more complicated than those presented herein. 


THEORY OF METHOD 

The method of analysis presented in this section determines the 
velocity distribution along the profiles of blade elements on surfaces 
of revolution. 
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Preliminary Considerations 

Assumed nature of flow . - In this section certain preliminary 
assumptions are made concerning the three-dimensional flow of an ideal 
compressible fluid through an arbitrary impeller passage between blades 
such as shown in figure 1. In general, the fluid is free to follow 
whatever path the pressure and inertia forces require of it. If, 
however, it is assumed that the number of blades in the impeller 
approaches infinity, the space between blades approaches zero and the 
path of the fluid is restricted to the curved, mean surface of the 
blade. (The blades become very thin so that the two surfaces of each 
blade approach a mean surface.) Under this assumption of axial 
symmetry the fluid motion is reduced from a general three-dimensional 
motion to a two-dimensional motion on the curved, mean blade surface. 
The streamlines of this two-dimensional motion can be projected on the 
meridional (axial -radial) plane, as shown in figure 2. Ruden (refer- 
ence 10) has shown that, provided that the blades are not too widely 
spaced, axial-symmetry solutions give a good picture of the mean flow 
between blades. 

For finite blade spacing, flow conditions vary between blades in 
the circumferential direction about the axis of the impeller. In order 
to investigate this blade-to-blade variation, it is assumed that the 
motion of any fluid particle bounded by adjacent streamlines in the 
meridional plane (fig. 2) is restricted to the annulus generated by 
rotating these adjacent streamlines about the axis of the impeller. 

If the adjacent streamlines are sufficiently close together, flow 
■conditions in the annulus can be considered uniform normal to a mean 
surface of revolution in the annulus. Thus the fluid motion is reduced 
to two-dimensional flow on the mean surface of revolution (fig. 3) 
generated by rotating the center line between the adjacent streamlines 
in the meridional plane (fig. 2) about the axis of the impeller. 

Blade-to-blade solutions of this type may be obtained for every 
mean surface of revolution generated by the center lines between 
adjacent streamlines in the meridional plane. Therefore, flow condi- 
tions can be determined throughout the passage between blades. The 
resulting quasi three-dimensional solution is obtained by the combina- 
tion of two types of two-dimensional solution, axial-symmetry solu- 
tions in the meridional plane and blade-to-blade solutions on surfaces 
of revolution. Such a combination of solutions prohibits the possi- 
bility of a corkscrew path, which the fluid might follow in an exact 
three-dimensional solution, but it can be expected to give a better 
picture of the flow than does any two-dimensional solution alone. 

The method of analysis just described is accomplished in two 
phases, axial- symmetry solution and blade-to-blade solutions. Only 
the second phase, blade-to-blade solutions, will be considered in this 
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report. The shape and the distribution of meridional streamlines in the 
axial -radial plane are assumed to be known from an axial -symmetry- 
solution (reference 11, for example). Thus, for a hlade-to-hlade 
solution in the annulus generated about the axis of the impeller hy any 
two adjacent meridional streamlines (fig. 2), the shape of the mean 
surface of revolution (fig. 3) is known from tie shape of the center 
line between the adjacent meridional streamlines, and the variation in 
height of an elementary fluid particle (fig. 2) as it moves along the 
mean surface of revolution from the impeller inlet to the impeller tip 
is known from the variation in spacing of the adjacent meridional 
streamlines . 

Coordinates . - The cylindrical coordinates R, Q, and Z are 
shown in figure 3. (All symbols are defined in the appendix. ) These 
coordinates are dimensionless, the linear coordinates R and Z having 
been divided by the impeller-tip radius rp. The coordinate system is 
oriented with the Z-axis along the axis of the impeller. The coordi- 
nates are fixed relative to the impeller, which rotates with the 
angular velocity <n in the positive direction (right-hand rule) about 
the Z-axis, as shown in figure 3. 

An infinites imal distance dS in the direction* of flow (that is, 
coinciding with the velocity vector) has components dR, Rd0, and dZ 
(fig. 3). The projection of dS on the meridional plane is given by 
dM in figure 3. The infinitesimal distances dS and dM help to 
define two angles a and £ where, from figure 3, 


dR = dM sin a 

(la) 

dZ = dM cos a 

(lb) 


and 


dM = dS cos p 

(2a) 

Rd0 = dS sin 0 

(2b) 


The angle a (fig. 3) is determined by tangents to the center line, 
between adjacent meridional streamlines, that generates the surface of 
revolution. The angle (3 (fig. 3) is the flow direction on the sur- 
face of revolution measured from a meridional line. From equation (2a) 



because dS and dM are always positive and finite. From equa- 
tions (la) and (lb), 


0< a — f 
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because for impellers of centrifugal compressors dR and dZ will be 
considered positive (or zero). 

Fluid strip . - A fluid strip of infinitesimal width dM lies on 
the surface of revolution and extends across the passage between blades 
along a line of constant R. A developed view of the fluid strip is 
shown in figure 4. The fluid strip has dimensions dM and RA0 where 
the angular width of passage between blades A 9 is defined by 

A0 = 0 t -0d (3) 

in which the subscripts d and t refer to the driving and trailing 
faces of the blades, respectively (left and right walls of the channel 
between blades in fig. 4). The height ratio H of the fluid strip is 
defined as the ratio of the incremental height Ah (fig. 2) of the 
fluid strip at radius R to the incremental height (Ah)>p of the 
fluid strip at R = 1.0. This height ratio is completely determined 
along a mean surface of revolution by the spacing between the adjacent 
streamlines in the meridional plane (fig. 2) . 

Velocity components . - The relative velocity Q on a surface of 
revolution has components and Q q in the dM and d0 directions, 

respectively, (fig. 3). These velocities are dimensionless, having 
been divided by the absolute stagnation speed of sound c Q upstream of 
impeller, where 


c 0 2 = rgi®o ( 4 ) 

in which R is the gas constant, y is the ratio of specific heats, 

T is the static (stream) temperature and where the subscript o refers 
to stagnation conditions upstream of the impeller. The tip speed of 
the impeller is likewise dimensionless and equal to the impeller-tip 
Mach number M^, which is defined by 

(JOCTm 

Mr = - — (5) 

c o 

Thus, the tangential velocity of the impeller at any radius R is 
equal to RMip and the absolute tangential velocity of the fluid is 
equal to (RM^ + Qg). From figure 3 

% = Q cos 0 (6) 

and 


ft. 


Qg = Q sin 0 


(V) 



6 


NA2A TN 2421 


4 


Thermodynamic relations . - From the general energy equation and 
from the isentropic relation between temperature and density, the 
density ratio p/p Q is related to the relative velocity Q by 



(RM t ) 2 -Q 2 -2M t X u | 


1 

r-i 


(Q) 


where the subscript U refers to conditions upstream of the impeller 
and where A is the whirl ratio (absolute moment of momentum divided 
by r T c 0 ) given by 


~K = R(RMrp + Qq) 


( 9 ) 


Development of Method 

Assumptions . - Before outlining the method of analysis it is con- 
venient to discuss the major assumptions. Consider the fluid strip 
in figure 4. Along the infinitesimal distances bounding the fluid 
strip at the driving and trailing blade surfaces, the velocities may be 
considered constant and equal to Q 4 and Q-^, respectively, and the 
flow directions may be considered constant and equal to P^ and p^. 
Along the lines of constant R bounding the fluid strip in figure 4, 
the velocity varies in some unlmown manner from to Q-^ and the 

flow direction varies from P^ to P^. In this report it is assumed 
that the average values of Q and P along lines of constant R may 
be used to satisfy the conditions of continuity and absolute irrota- 
tional motion. The average value of Q is assumed to be given by 


Q. 


_ Qd + Qt 


-av 


( 10 ) 


and, for the average value of P is assumed to be given by 


3d + 3t 


■'av 


(id 


Also, for R X <R< 1 . 0 , 


sin P av = A + BR + CR 2 


( 12 ) 
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where A, B, and C axe coefficients to he determined and where R x 
is the largest radius at which the fluid is considered to he perfectly 
guided hy the blades; that is, the radius at which the simplified 
analyses given in references 6 and 7 break down. From figure 10 of 
reference 7 the value of R x for A0 equal to rt/lO and sin a 
equal to 1.0 is about 0.8. For other values of Ad and sin a, the 
value of Rjj. can be estimated from 

In R x _ in 0.8 _ 

(A0 ) sin a rt/10 

where a is the average value over the interval R x — R 5. 1.0. Equa- 
tion (13) is based upon an extension of the work in reference 1 where 
for impellers with straight blades it is shown that the flow conditions 
in one impeller can be correlated with the flow conditions in another 

**| 

impeller at the same value of > -rrrr : . In reference 1, the 

(A0) sin a ’ 

impeller-tip Mach number and the compressor flow rate were found to have 
a negligible effect on the value of R x . 

Outline of theory . - Fluid strips such as shown in figure 4 exist 
at all radii along the surface of revolution. From the assumptions of 
this analysis there are three unknowns (Q^, (fy., and ) for each 

fluid strip. These unknowns can be determined by the simultaneous 
solution of equation (10) and the equations of continuity and zero 
absolute circulation for flow across the fluid strip. Equations for 
the distribution of velocity along the blade profile on a surface of 
revolution will be developed in this report. 

Zero absolute circulation . - In the absence of entropy gradients, 
which result from shock, viscous dissipation, heat transfer, and so 
forth, the absolute circulation around the fluid strip in figure 4 is 
zero so that 

(Q + RM T sin P) d - (Q + RMt sin P) t + 



_d^ 

dM 


(RMt + Qe) av BAB 


dM = 0 


(14) 


where (RMj + Qg) av is the average absolute tangential velocity and 
where from trigonometric considerations of the velocity triangles 
(fig. 5) Q + RMj sin (3 is the absolute velocity component along the 
blade surface. From equations (7) and (14) and from the assumptions 
that Q and p equal Q^ v and P av , respectively, in the passage 
between blades, 
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Qd 


Qt 


+ RMrj (tan P d - tan P-^) + 


cos Pa cos P^ 

A ^(rMt + S in p ay ) RA©] = 0 


( 15 ) 


Finally, from equations (10) and (15) 


cos Pa cos p t I 2 Qay 

Qt = ft, 7 Vno R n i nng~~RT + ( ta “ P d - ^ P t ) + 


COS fi-k + COS I COS 

(RMr + Qav sin p av ) R(Ae)]' 


d 

dM 


and from equation (ll) 


Qd ~ ^Qqv “ Qt 


(16) 


(17) 


If Qav and p av are known, Qa and can be determined from 

equations (16) and (l7). 

Average velocity Qav- - From continuity considerations of the flow 
across the fluid strip in figure 4, 

p av 


Aw = P 


° P, 


( c o^av ^ cos p av (^T H ( r T^^ ^ 


from which 


Qav - 


J£_ 


Pav 


cos p av HR 


A9 

IaqTt 


where the flow coefficient <p is defined by 


cp = 


1 Aw 
p o c o (^ & )t 


(18) 


(19) 


in which Aw is the incremental flow rate through the passage between 
two blades on the surface of revolution and (Aa)rp is the incremental 
flow area (between two blades) normal to the direction of at the 

impeller tip 


( Aa ) qi — r»p ( A0 ) rp (Ah ) ip 


(19a) 


J/f/ 
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The flow rate per unit flow area at the impeller tip Aw/ (Aa ) ^ is 
known so that the flow coefficient cp can be determined by equa- 
tion (19). The density ratio P av /p 0 is given by equation (s) with 

Q equal to Q av 



Therefore, the velocity Qgy can be determined by the simultaneous 
solution of equations (18) and (20) provided that the average flow 
direction P av is known. 

Average flow direction . - In the passage between blades the 
average flow direction is assumed equal to the average blade direction 
(equation (ll)) except near the blade tip (R x < R < 1.0) where 
sin 3 av is given by equation (12) . The exact variation in 8 av 
with R in the interval R x < R <1.0 could be represented by an 
infinite series. However, because the variation in sin P av with R 
will not, in general, contain an inflection point, a parabolic varia- 
tion in sin 3 av with R has been assumed and only the first three 
terms of the infinite series retained. The constants A, B, and C 
in equation (12) are determined from: 

(1) (sin 3 av ) x = A + BR X + CR X 2 



and 

(3) (sin 3 av )ij — A + B + C 

so that 
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where (sin p av ) x n — — ^ are known from the blade geometry 

at R x , and where (sin P av )p ls determined from the slip factor p, 
which is defined by (reference 6) 


H = 1 + 


(Qav s ^ n ^av^T 

Mip 


so that 


(sin P av ) T = 


Mt (p-l) 
(Oav^T 


( 22 ) 


The slip factor p is assumed to be known, or can be estimated, as a 
result of the work presented in references 1 and 7, for example. 
(Further discussion on R x and the slip factor p is given later in 
this report.) The velocity (Qav)T equation (22) is obtained from 
equations (18) and (20) with R and H equal to 1.0 and with 

(cos P av )p replaced by ^ 1 - (sin P av )p 2 where (sin p av ) T is given 
by equation (22): 



Equation (23) is solved for (Og^T ty trial and error. Therefore, 

is determined as a function of R (or M) by equations (ll), (12) , 

cl V 

and (21). The velocities Q<p, and Qgy are determined as func- 

tions of M (or R) from equations (16), (17), and (18). (The last 
term of equation (16) is determined from the slope of 
(RMqi + Qg V sin P av ) R (A0) plotted against M. ) 


APPLICATION OF METHOD 

The following outline of the numerical procedure is given for the 
general case of a mixed-flow impeller with arbitrary hub and shroud 
contours in the meridional plane (fig. 2) and arbitrary blade shape 
(curvature and thickness distribution) on surfaces of revolution. It 
is assumed that the surfaces of revolution are known, that is have been 
generated by the center lines between adjacent meridional streamlines 
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obtained from an axial symmetry solution (reference 11, for example). 

The following outline of the numerical procedure refers to any one of 
these surfaces of revolution. 

Specified conditions . - The following quantities are specified: 

(1) Flow coefficient cp (defined by equation (19) in which 

p , c Q , Aw, and (Aa)^ are known quantities) 

(2) Impeller-tip Mach number (defined by equation (5)) 

(3) Whirl ratio X^ upstream of impeller (defined by equation (9)) 

(4) Ratio of specific heats y 

(5) From the shape of the center line between adjacent meridional 

streamlines that generate the surface of revolution, 

R = R(M) 

and 

a = a(M) 

where the distance M along a meridional line on the sur- 
face of revolution is arbitrarily equal to zero at the impel- 
ler tip and decreases toward the impeller inlet 

(6) From the spacing of the adjacent meridional streamlines, 

H = H(M) 

Variation in sin P av . - The variation in sin g av with R near 
the impeller tip (R x < R < 1.0) is determined as follows: 

(l) Compute the value of R x by equation (13) . If sin a varies 
in the region H x < R <1.0, as it generally does, the 
average value of a in this region is used in equation (13), 
and because the average value of a varies with the value 
of R x , equation (13) must be solved by trial and error. 
However, because the value of a does not generally vary 
greatly in the region R x < R < 1.0, a satisfactory value 
of R x could be obtained from equation ^13) using the 
average value of a obtained from an initially assumed 
value of R x . Also, equation (13) was developed from 
information (references 1 and 7) relating to blades that 
are not designed to unload at the tip. If the blades 
being considered were designed to unload at the tip, the 
direction of the mean flow path near the impeller tip 
would deviate less from the mean blade direction and the 
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value of Rx would be somewhat greater than that given by 
equation (13) . The value of R x is not especially criti- 
cal and in these cases, with sufficient experience, it can 
probably be estimated accurately enough from the radius at 
which an assumed path of the mean streamline (sketched by 
experience) deviates appreciably from the mean direction 
of the blade profile on the surface of revolution. 

(2) Est ima te the value of p, or obtain values from references 1 

and 7. The values of p given in references 1 and 7 were 
obtained for blades that are not designed to unload at the 
tip. If the blades being considered were designed to unload 
at the tip, the direction of the mean flow path at the impel' 
ler tip would deviate less from the mean blade direction 
at the tip and the value of p would be somewhat greater 
than that indicated in references 1 and 7; that is, OWt 
would be more nearly equal to 77 (P^ + P-t) 1 !* The value 

of p, like the value of R x , is not especially critical 
and in these cases it can probably be estimated accurately 
enough from the assumed shape of a mean stre aml ine (sketched 
from experience) between blade-element profiles on the sur- 
face of revolution. 

(3) Compute (Qav)T from equation (23) by trial and error. 

(4) Compute (sin P av )i from equation (22). 

(5) Compute coefficients A, B, and C from equations (2l). 

(6) Compute sin p over the interval R,< R < 1.0 by equa- 

tion (12). 

The variation in p av for R less than R^. is given by equa- 
tion (ll). This equation is assumed to be valid downstream to the 
impeller inlet. If the angle of attack at the impeller inlet is zero, 
the assumption is probably good. If the angle of attack is small, the 
error involved is probably small and could be partly corrected by 
est ima ting the path of the mean streamline between blades in this 
region. For large angles of attack, the stagnation point on the blade 
surface may exist well inside the impeller passage and the ideal flow 
is reversed along the blade surface downstream of this point. Under 
these conditions, the method of analysis does not apply near the 
impeller inlet but because of the high blade solidity it does apply 
elsewhere in the impeller. 
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Average velocity Q av . - The average velocity Qay at each value 
of M (or R) is determined "by equation (18) in which p av /p Q is given 
hy equation (20). Because P av /p 0 also contains Q av , the simultan- 
eous solution of equations (13 ) and (20) must he hy trial and error. 

A suggested procedure is first to compute P av /p Q assuming that Q av 
in equation (20) is zero. This value of P av /p 0 is then used to com- 
pute Q ay hy equation (18) . The process is repeated each time using 
the new value of to compute P av / P Q until the value of Q av 

converges . 

Velocities on hlade surfaces, and - The velocities on the 
hlade surfaces at each value of M (or R) are determined hy equa- 
tions (16) and (17). The last term in equation (16) is obtained from 
the slope of (EM^ + Q^y sin 3 av ) R (A0) plotted against M. 

Finally, the static (stream) pressure p corresponding to the 
relative velocity Q at any radius R is given hy 

£ - (fj " j 1 + [(i%) 2 -S 2 -2M T \,]j ?T (24) 


NUMERICAL EXAMPLES 

The approximate analysis method developed in this report is 
applied to eight examples for which relaxation solutions of the exact 
differential equation for two-dimensional compressible flow in impel- 
lers of centrifugal compressors are given in references 1 and 7. 
Although these examples are for radial- or conic-flow surfaces and not 
for arbitrary surfaces of revolution, they cover a fairly wide range 
of design and operating variables so that a comparison of the veloci- 
ties (on the hlade surfaces) obtained hy the relaxation solutions and 
hy the approximate analysis method should serve as a check on the 
validity of the approx im ate method. 

Types of impeller . - The eight numerical examples are for two- 
dimensional radial -flow inn =>Uers for which a is equal to 90° and the 
surfaces of revolution are flat planes normal to the axes of the 
impellers. The impellers (fig. 6) contain a finite number of thin 
straight (p^ = p-^ = 0) or logarithmic-spiral (p^ = p^ = constant) 

blades, and the flow area normal to the direction of Qm is constant 
so that HR equals 1.0. Only the critical flow region toward the tip 
of the impeller was investigated (0.70< R < 1.0). The diffuser vanes 
(if any) and the inducer vanes were assumed to he far enough removed 
not to affect the flow in this region. 
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Design and operating variables . - The following design and opera- 
ting variables were specified for the eight examples: 


Example 

cp 

% 

A0 

tan p 

Fluid 

g 

(a) 

0.5 

1.5 

2rt/30 

0 

Compressible 

0.934 

(*> 

.7 

1.5 

2rt/ 30 

0 

Compressible 

.937 

(c) 

.9 

1.5 

2rt/30 

0 

Compressible 

.938 

(4) 

.5 

2.0 

2jt/30 

0 

Compressible 

.935 

(e) 

.5 

1.5 

2rt/ 20 

0 

Compressible 

.899 

(f) 

.5 

1.5 

2rt/ 20 

-0.5 

Compressible 

.834 

(g) 

.5 

1.5 

2*/20 

-1.0 

Compressible 

.768 

(H) 

.5 

1.5 

2jt/20 

0 

Inc ompr e s s ibl e 

.892 


The whirl ratio upstream of the impeller Ay was zero and for the com- 
pressible fluid the ratio of specific heats y was 1.4. The value 
of the slip factor p. given in the table was obtained from the relaxa- 
tion solutions and was also used to compute (sin p av )gi i* 1 equa- 
tion (22). For the incompressible example, the speed of sound c D con- 
tained in the definitions of Q, <p, and Mg’ is a fictitious quantity 
(constant) considered equal to the upstream stagnation speed of sound 
of the compressible-flow examples with which the incompressible -flow 
example is compared. 

Results . - The results of the comparison between the relaxation 
solutions and the approximate method of analysis are shown in figure 7 
for the eight examples. The velocities and fy. are plotted 

against R for the relaxation solutions and for the approximate solu- 
tions. The average velocity Qg V used in the approximate method of 
analysis to obtain Q<g and is also plotted. 

The effect of increasing the flow rate (flow coefficient) on the 
agreement between the relaxation solution and the approximate solution 
is shown in figures 7(a) to 7(c). The agreement appears equally good 
for all flow rates. In view of the relative simplicity of the approxi- 
mate method of analysis, the agreement is considered entirely satis- 
factory ' everywhere except in the immediate vicinity of the impeller tip 
where the results of the approximate method of analysis do not follow 
the rapid unloading of the blades. This rapid unloading is characteris- 
tic of blades that are not designed to unload at the tip. If the 
blades were designed to unload, the agreement between relaxation solu- 
tions and approximate solutions should be better. In any event the 
disagreement is serious only over the last 2 percent of impeller -tip 
radius. The failure to unload at the impeller tip will be observed in 
most of the remaining solutions. 
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The following figure comparisons indicate: figures 7(a) and 7(d), 

the effect of increasing impeller -tip Mach number; figures 7(a) 
and 7(e), the effect of increasing angular width of passage between 
blades (A0); figures 7(e), 7(f), and 7(g), the effect of larger 
negative blade angles P; and figures 7(e) and 7(h), the effect of 
compressibility. In figure 7(h) , the peculiar humps in the velocity 
distribution obtained by the approximate method of analysis indicates 
that for incompressible flow the blades start to unload at a lower 
value for than that given by equation (13) . 

In view of the relative simplicity of this approximate, but rapid, 
method of analysis, the agreement between the relaxation solutions and 
the approximate solutions is considered good in all cases investigated; 
that is, over fairly wide ranges of flow rate, impeller -tip speed, 
blade curvature, and number of blades. 


SUMMARY OF RESULTS 

A rapid approximate method of analysis was developed for deter- 
mining the velocity distribution on impeller blades of centrifugal 
compressors. In eight numerical examples the velocities obtained by 
the approximate method of analysis were compared with the more nearly 
correct values obtained by relaxation methods. In all cases, that is, 
over a fairly wide range of flow rate, impeller -tip speed, blade 
curvature, and number of blades, the agreement between velocities 
obtained by the approximate method of analysis and by relaxation 
methods was considered good. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 27, 1951. 
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APPENDIX - SYMBOLS 

The following symbols are used in this report: 

A,B,C coefficients defined by equation (2l) 

c Q stagnation speed of sound upstream of impeller, equation (4) 

g acceleration due to gravity 

H height ratio of fluid strip normal to surface of revolution, 

Ah/(Ah) T 

M distance along meridional line on surface of revolution 

(dimensionless, expressed as ratio of impeller tip radius rrp) 
(fig. 3) 

impeller tip Mach number, equation (5) 

p static (stream) pressure 

Q relative velocity on surface of revolution (dimensionless, 

expressed in units of the stagnation speed of sound upstream 
of impeller c Q ) (fig. 3) 

R cylindrical coordinate (dimensionless, expressed as ratio of 

impeller-tip radius rtf) (fig. 3) 

r T impeller-tip radius 

S distance along streamline on surface of revolution (dimension- 

less, expressed as ratio of impeller-tip radius rip) (fig. 3) 

T static (stream) temperature 

Z cylindrical coordinate (dimems ionless, expressed as ratio of 

impeller-tip radius rip) (fig. 3) 

a slope of surface of revolution in direction of %, equa- 

tions (la) and (lb) (fig. 3) 

p flow direction on surface of revolution, equations (2a) and (2b), 

(fig. 3) 

X ratio of specific heats 

(Aa)gi incremental flow area between two blades and normal to the 
direction of at impeller tip, equation (l9a) 
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Ah incremental height of fluid strip on surface of revolution 

Aw incremental flow rate between two blades on surface of revolu- 

tion 

A 9 angular width of passage between blades, radians unless other- 

wise specified, equation (3) 

9 cylindrical coordinate, radians unless otherwise specified, 

(positive about Z-axis according to right-hand rule) (fig. 3) 

A whirl ratio, equation (9) 

p slip factor, equation (22) 

p static (stream) weight density of fluid 

cp flow coefficient, equation (19) 

cn singular velocity of impeller (in direction of positive 9 ) 

Subscripts: 

abs component of absolute velocity along blade surface 
av average 

d driving face of blade (blade surface in direction of rotation) 

(fig. 4) 

M component along meridional line on surface of revolution 

o absolute stagnation condition upstream of impeller 

R,0,Z components in positive R-, 0-, Z-directions, respectively 

T impeller tip 

t trailing face of blade (blade surface opposed to direction of 

rotation) (fig. 4) 

U upstream of impeller 

position along meridional line on surface of revolution at which 
the assumption of perfect guiding of fluid by blades is con- 
sidered to break down 


x 
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Figure 2. - Streamlines in meridional plane for axial -symmetry solution 
of flow through impeller of figure 1. 
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(b) Example (b): flow coefficient 0.7; other parameters same as example (a). 

Figure 7. - Continued. Variation in velocity along blade surfaces as obtained by relaxation 
methods (references 1 and 7) and by approximate method. 
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(d) Example (d): impeller-tip Mach number M T , 2.0; other parameters same as example (a). 

Figure 7. - Continued. Variation in velocity along blade surfaces as obtained by relaxation 
methods (references 1 and 7) and by approximate method. 
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-Relaxation solution 
-Approximate method 




(e) Example (e): angular width of passage A 0, 18 ; other parameters same as example (a). 

Figure 7 . - Continued. Variation in velocity along blade surfaces as obtained by relaxation 
methods (references 1 and 7) and by approximate method. 





(f) Example (f): blade angle fj, tan” 1 (-0.5); other parameters same as example (e). 

Figure 7 . - Continued. Variation In velocity along blade surfaces as obtained by relaxa 
tion methods (references 1 and 7) and by approximate method. 
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(g) Example (g) : blade angle tan“ 1 (-1.0); other parameters same as example (e). 

Figure 7. - Continued. Variation in velocity along blade surfaces as obtained by relaxa 
tion methods (references 1 and 7) and by approximate method. 












